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Abstract 

Introductory talk at the Session on Correlations and Fluctuations 
of the XXVIIth International Symposium on Multiparticle Dynamics 
held at Delphi, Greece from 6th till 11th of September 1998. 

1 Introduction 

The purpose of this talk, as I understand it, is to introduce briefly the subject 
of the Session to the participants who are not experts in the field. In trying 
to do this I shall heavily borrow from my recent summary of the Matrahaza 
workshop |I| . Since I am now much less constrained by the program of the 
meeting, however, this account shall reflect more adequately my personal 
views on the subject. I restrict myself to the four issues: 

(i) Bose-Einstein interference; 

(ii) Intermittency; 

(iii) QCD and multiparticle correlations; 

(iv) Event-by-Event fluctuations. 
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2 Bose-Einstein interference 



The discussion of correlations in multiparticle production is at present largely 
dominated by effects of the Bose-Einstein interference. Let me thus start by 
a brief reminder what is this all aboutQ. 

The practical problem we face can be formulated as follows: given a 
calculation (or a model) which ignores identity of particles, the question 
is how to "correct" it in order to take into account the effects of quantum 
interference (which is the consequence of identity^). Let us thus suppose that 
we have an amplitude for production of N particles M$\q) (q = qi, ...gjv), 
calculated with the identity of particles being ignored. The rules of quantum 
mechanics tell us that, to take the identity of particles into account, we 
have to replace M$\q) by a new amplitude M^{q) which is a sum over all 
permutations of the momenta (qi, ...gjv) 

M$\q)-,M N (q)=Y, M NM- (1) 

p 

This would be the end of the story if particle production was described 
by a single matrix element. In general, however, we have to average over 
parameters which are not measured and therefore the correct description of 
the multiparticle final state is achieved in terms of the density matrix 

pSW) =^M§\q,u)M!^(M, (2) 

rather than in terms of a single production amplitude. The sum in (Q) runs 
over all quantum numbers u which are not measured in a given situation. 

(q, q') gives all available information about the system in question. At 
this point it is useful to note that, when tranformed into (mathematically 
equivalent) Wigner representation 

W N (q,x) = J d(Aq)e^p^(q,Aq) (3) 
lr The physics of Bose-Einstein correlations was recently extensively reviewed by G.Baym 

| 

2 It should be understood that this problem is very common in quantum mechanical 
calculations, as illustrated, e.g., by evaluation of Feynman diagrams. I would like to thank 
J.Pisut and K.Zalewski for discussions of this question. 
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(q — {q + Q ,/ )/2; = q — q') it gives information about the distribution of 
momenta and positions of the particles (see, e.g., for a discussion of this 
point). 

Using (|l]) and (0) one easily arrives at the formula for the corrected (i.e., 
with identity of particles taken into account) density matrix P7v(g, g') and 
one finally obtains the observed multiparticle density 

^) = ^Ep ( °Wp') (4) 
iV • p,p' 

where the sum runs over all permutations P and P' of the momenta (gi, ...gjv). 
The factor appears because the phase space for N identical particles is 
iV! times smaller than the phase space for iV non-identical particles. The 
formula (f|) is in common use[] and is the basis of our further discussion. 

2.1 A theoretical laboratory: independent particle pro- 
duction 

The case of independent particle production is an attractive theoretical lab- 
oratory which, although not expected to describe all details of the data, 
reveals -nevertheless- some generic features of the problem. This was first 
recognized by Pratt J|. In terms of the density matrix, the independent 
production means that the density matrix factorizes into a product of single- 
particle density matrices 

pSW) = P { °\Qi,q'x)p {0) (q2,q' 2 )....p (0) (qN,q' N ) (5) 
and that the multiplicity distribution is the Poisson one 

P<°)(JV) = e -*— . (6) 

It turns out ||, []] that in the case of a Gaussian density matrix the 
problem can be solved analytically. The main results (valid also in the general 
case of an arbitrary density matrix 0) can be listed as follows. 

3 Using the symmetry properties of the density matrix, the double sum in (^) can be 
reduced to a single sum. The factor -™ is then absent. 
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(a) All correlation functions K p (qi, q p ) and the single particle distri- 
bution Q(q) can be expressed in terms of one (hermitian) function L(q, q') = 
L*(q', q) of two momenta: 

Q(q) = L(q, q); K 2 (q 1: q 2 ) = L(q h q 2 )L(q 2l q%)] 
K 3 (qi, 92, 93 ) = L(q x , q 2 )L(q 2 , q 3 )L(q 3 , qi) + L(q 1 , q 3 )L(q 3 , q 2 )L(q 2 , q x ), (7) 

and analogous formulae for higher correlation functions. 

(b) At very large phase-space density of particles, the distribution ap- 
proaches a singular point representing the phenomenon of Bose-Einstein 
condensation: almost all particles populate the eigenstate of p^\q.,q') cor- 
responding to the largest eigenvalue. The resulting multiplicity distribution 
is very broad (almost flat) so that, e.g., the probability of an event with no 
single 7T° produced is non-negligible 0. 

It should be not surprizing if the very restrictive condition of independent 
production, as expressed by (|5]J|) is not realized in nature. Nevertheless the 
comparison of the relations (|?|) with data is interesting, since they are a sort 
of reference point allowing to judge if the observed multiparticle correlations 
are "large" or "small" with respect to the observed two-body correlations. 
The existing evidence leads to rather interesting, although controversial, con- 
clusions. At the Matrahaza meeting, Lorstad [TI| demonstrated that there 



are practically no genuine three-particle correlations^ in S-Pb collisions at 
CERN SPS. Since the two-particle correlations are clearly visible, this obser- 
vation is not easy to reconcile with Eq.(7). It was earlier shown by Eggers et 
al [12] that the UAl data are also in contradiction with (0), although in this 



case the 3-body correlations seem to be too large to satisfy (|7|). On the other 
hand, it was shown recently by Arbex et al that the NA22 data agree 



well with (|7|). This striking difference between the behaviour of heavy ion 
and "elementary" collisions is certainly very interesting and deserves further 
attention. 

We cannot thus consider the results obtained from (|5],^D to be a realistic 
description of the data. Nevertheless, the main conclusion about the possi- 



4 This effect was also considered in connection of the possible production of the Dis- 
oriented Chiral Condensate ||. The present argument adds another obstacle on the 
difficult road to observation of DCC. 

5 The importance of the absence of 3-particle correlations in heavy ion collisions was 



emphasized already some time ago |11 
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bility of Bose-Einstein condensation remains an interesting option which is 
worth a serious consideration [1T41 |15 . 



2.2 Monte Carlo simulations 

In this situation, the practical method to study the effects of BE symmetriza- 
tion on particle spectra is to implement it into the Monte Carlo codes. A 
"minimal" method of performing this task was suggested some time ago [[|. 
The idea is to take an existing code (which reproduces the distribution of 
particle momenta, i.e. the diagonal elements of the density matrix) and to 
introduce an ansatz for the off-diagonal elements of the multiparticle density 
matrix ([?])[]. Each event generated by the MC code is then given a weight 
which is calculated as the ratio of symmetrized distribution [Eq.(f|)], and the 
unsymmetrized one. In this way the modification of the original spectra is 
kept at the minimum. 

A practical realization of this idea (in its simplest version) has been de- 
velopped by the Cracow group |16| and shall be presented by Fialkowski at 



this meeting. They propose the unsymmetrized density matrix in the form 

N 



P %\q,q') = PN{qmw{qi-ti (8) 

i=i 



where Pn{q) is the probability of a given configuration obtained in JET- 
SET and if is a Gaussian. This prescription does not modify the diagonal 
elements of the unsymmetrized density matrix {w(0) = 1) and, moreover, 
does not introduce any new correlations between emission points of the pro- 
duced particles (when transformed into Wigner representation, Eq.([|), the 
product Y[ w (<li ~ q'i) becomes the product Y[w(xi)). Thus (H) can indeed 
be considered as a minimal modification of the existing code. The autors 
find that this prescription represents well the existing data on two-particle 
correlations and that they can recover the experimental multiplicity distri- 
bution by a simple rescaling with the formula P{N) — > P(N)cV N , without 
the necessity of refitting the JETSET parameters. 

They also studied production of pair of W bosons at LEP II |L6| and 
found fairly strong effects of quantum interferrence. One may note, however, 



6 As seen from this corresponds to introducing an - a priori arbitrary - distribution 
of particle emision points in configuration space. 
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that in present version of the model the position of particle emission point 
is not correlated with its momentum, whereas such correlation is likely to be 
present in reality. Consequently, the obtained results may be overestimating 
the effect. This deficciency is easy to repair || but on the prize of introducing 
more parameters. 

A more fundamental approach has been pursued since some time by An- 



dersson and Ringner [17] and shall be presented here by Todorova-Nova. It 



is based on the paper by Andersson and Hoffman |Tj|. In this case the "un 



corrected" matrix element represents the decay of one Lund string which 
is then symmetrized according to the procedure explained in introduction. 
Two particle correlations are well described and several interesting effects are 
predicted. Among them: (a) the longitudinal and transverse correlations are 
expected to be different because they are controlled by two different physical 
mechanisms; (b) Three particle correlations are predicted non- vanishing and 
were actually calculated; (c) WW production was studied and no significant 
mass shift is expected; (d) No multiplicity shift in the W decay is predicted. 

This last conclusion is a consequence of the fact that, in case of more 
than one string present in the final state, no symmetrization between parti- 
cles stemming from different strings is performed. This corresponds to the 
assumption that the strings are created at a very large distance from each 
other. One thus may expect that in a more realistic treatment some multi- 
plicity shift should be present [] 

Finally, let me add that in both [16| and |17| the "interconnection effect" 
(which has tendency to reduce the multiplicity) is neglected. The full 
phenomenological analysis of the data is therefore certainly more compli- 
cated, as we shall hear from de Jong. 



2.3 Probing the space-time structure 

Much attention is also devoted to the information one may obtain from the 
data on quantum interference about the space-time structure of the multi- 
particle system created in the collision [c.f.Eq.(fJ)]. Although such analyses 
have a somewhat limited scope, as they (i) provide only information about 
the system at the freeze-out and (ii) require several additional assumptions 
- they give nevertheless a unique opportunity to investigate this problem. 

7 A contribution to this problem was presented recently by B.Buschbeck et al. [|l9| 
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Most of the caveats are thus usually postponed to the future (and better 
data) and the analysis is carried on. 

The recently presented investigations were based on the hydrodynamic 
approach. Some of them were discussed here already during the Session 
on Heavy Ion Interactions. The main features are: (i) The shape of the 
particle emission region is consistent with the in-out scenario of Bjorken |H] 



(ii) The longitudinal size of the "fireball" from which a bulk of particles 
are emitted is several times larger in heavy ion collisions than in hadron- 
hadron interactions; (iii) Particle emission process starts rather late in heavy 
ion collisions (after about 4 fm in S-Pb interactions []), as compared to the 
elementary collision where it happens immediately after collision []; (iv) the 
emission process, once started, does not last for a long time: less than 2 fm 
for elementary and about 3 fm for S-Pb interactions. 

These features clearly indicate that a heavy ion collision is indeed followed 
by creation of an longitudinally expansing "fireball" in which some kind of 
matter is "boiling" for a considerable time. Once it is sufficiently cooled, 
however, its decay is rather fast. It is presumably too early to claim that 
it is formed from the quark-gluon plasma but nevertheless this behaviour is 
rather suggestive. 



3 Intermittency 

Intermittency |j22j| , postulates scaling of the multiparticle spectra. A rather 
complete review of the subject is now available |2^| , therefore I shall restrict 
myself to few remarks expressing my personal view on the progress achieved 
in the last decade. 

(i) The scaling hypothesis can be formulated in many ways and, indeed, 
much work was devoted to improvements and generalizations of the original 
proposal [^2J expressed in terms of the (normalized) factorial moments 

F q =< n(n - l)...(n - q + 1) > / < n > 9 ~ (bin sizeY fq . (9) 

The result was an impressive progress in the developpment of more sophis- 
ticated tools which are much better suited for investigation of many, some- 
times very detailed, aspects of the problem. Let me particularly emphasize 

8 A more detailed description is given in JlJ . 
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the importance of the correlation integrals, first introduced in Their use 
was decisive in improving accuracy of the data and thus to substantiate the 
evidence for the effect. 

(ii) On the experimental side, the analysis of high precision data (particu- 
larly those of NA22 and UAl experiments) allowed to establish -beyond a rea- 
sonable doubt- a close connection between intermittency and Bose-Einstein 
correlations, as suggested [^, Q almost immediately after first experimen- 
tal evidence for increasing factorial moments. This observation allowed to 
understand the scaling of the momentum spectra as a reflection of the -more 
fundamental- scaling in configuration space ^7|. I have impression that the 



importance of this fact is not yet fully appreciated. 

(iii) Recently, a general solution of the model of multiplicative cascade was 
obtained ||28|| . One can thus hope for a significant progress in understanding 



the scaling phenomenon. 

(iv) Finally, let me also mention another interesting development, namely 
the generalization of the notion of scaling by the idea of self-affinity p9| . I 
think it is an interesting direction to pursue and I hope that more data on 
the subject shall be soon available. We shall hear more about this from Liu. 

The major disappointment I see after all these years is that -in fact- no 
convincing theoretical basis was found for the phenomenon of intermittency, 



although it seems to be indeed an universal feature of particle spectra p3|. [30] 



[]. Does it mean that the simple effect one observes is a purely accidental 
result of summation of much more complex contributions? Perhaps. Never- 
theless, I am convinced that the search for a more fundamental reason of the 
apparent scaling in particle spectra is worth to continue. 



4 QCD and multiparticle correlations 

It is now rather well established that the average multiplicity and single 
particle spectra are well described by perturbative QCD supplemented with 
the principle of parton-hadron duality 



In my opinion, at present, the real challenge to the idea of parton-hadron 



duality is to explain the data on differential correlation functions | 27[| . Indeed, 



9 The second order phase transition was invoked by many authors pl[ | as a possible 
explanation. This is certainly a valid idea but it does not explain universality of the 
phenomenon. 
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it is hard to understand how the momenta of the produced hadrons can follow 
so closely the momenta of the created partons that the correlations between 
them are not washed outQ . Therefore a non-trivial extension of the principle 
of parton-hadron duality must be formulated in order to give quantitative 
meaning to perturbative calculations of multiple production. It was therefore 
rather recomforting to learn that indeed the predictions of perturbative QCD 



formulated some time ago p3| , are badly violated by the L3 data f34| . On 
the other hand, the same data are well described by the JETSET code. The 
conclusion is that the hadronization part is not correctly taken into account 
by the simple (naive?) parton-hadron duality. More about that later in this 
Session by Mandl. 

This is not to say that the subject is closed: Ochs [^] pointed out that the 
tested QCD calculations included several simplifying assumptions (the most 
important among them seems the neglect of energy-momentum consevation) 
and thus it is not obvious which part of the result is actually responsible 
for the failure. It is clear, nevertheless, that further work on these lines 
must seriously address the problem of parton-hadron duality and its range 
of application. 



5 Event-by-event analysis 

Event-by-event analysis clearly emerges as a next logical step in studies of 
multiparticle fluctuations. The subject is not yet well developped, however, 
and neither the physics nor methods sufficiently understood to define pre- 
cisely what we are really searching for. Therefore, I can only list a few ideas 
of potential interest. I am fully aware that some of them may not work and 
that others, more interesting, may well be proposed in the near future. 

There are two basic reasons why event-by event fluctuations attract at- 
tention. The first one, more spectacular, is to look for large deviations of 
some events from the average, with the hope of finding a hitherto unobserved 
effect. The second, more pragmatic, is to measure the distribution of a quan- 
tity defined for a single event and thus obtain additional information, helping 

10 This problem is much less serious if one considers only the integrated correlation 
functions. In this connection, see the discussion at the recent meeting on Correlations 
and Fluctuations, Matrahaza, June 1998 and contibutions to this session by Metzger and 
Chekanov. 
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to understand the physics of the process. This is well illustrated by multiplic- 
ity distribution which is the simplest event-by-event analysis one may think 
of. It was studied since long time Q and was of great help in understanding 
the physics of multiparticle production. 
Let me now go to my list: 

(i) Recently, Stodolsky |36| proposed to study fluctuations in transverse 
momentum (see also [[HJ). The idea is that, if the transverse momentum 
distribution in an event can be related to its "temperature", one obtains in 
this way the distribution of "temperatures" of the events. Now, if thermo- 
dynamics is a correct decription of the process in question, the fluctuations 



of temperature can be related to to heat capacity Cy of the system [[37 

= — (10) 

T 2 Cy 1 ' 

where k is the Boltzmann constant. This obviously may be very helpful in 
searching for phase transition. Even far from phase transition, however, a 
measurement of this kind can provide a lot of information on (a) the proper- 
ties of the system in question and (b) whether it is indeed close to thermo- 
dynamic equlibrium. To take a simple example: In case of ideal gas one has 
E = CyT, where E is the energy of the system. We thus obtain 

T 2 E K J 

The point is that both L.H.S. and R.H.S. of this equation can be measured 
and thus one may hope to estimate the deviation of the system from the ideal 
gas approximation. 

(ii) Another important issue was raised by Hwa He pointed out the 
essential difference between the determination of fractal parameters in case of 
dynamical systems and in case of systems of many particles. In the dynamical 
system one can generate the time sequence and thus estimate how fast the 
different trajectories diverge. In case of multiparticle systems we do not have 
a time sequence and thus we have to rely on patterns. The question in this 



11 Contributions related to multiplicity distributions shall be presented by Hegyi, Plosza- 
jczak and Blazek. 

12 Recently, the first results on temperature fluctuations were presented by NA49 coll. 
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case is: how different are the patterns of different events. Hwa proposed to 
measure the pattern of an event by the factorial moment associated with it. 
One can then ask the question how this measure fluctuates from event to 
event. Studying moments of this distribution provides a measure of event- 
to-event fluctuationF]. When they are considered as function of bin size, 
it is possible to define appropriate fractal dimensions which conveniently 
summarize the information. For the details the reader is referred to the 



original paper |39|. I personally feel that this is an important conceptual 
step in our thinking about the problem, although I am not fully convinced 
that the proposed measure cannot be improved. 

(iii) The studies of possible phase transitions in the multiparticle systems 
produced in high-energy collisions suggest that the fractal behaviour may 



strongly fluctuate from one event to another. It follows that it is essential 
to be able to study the fractal behaviour in event-by-event analysis. The 
feasibility of this program was investigated recently P5 |. It seems to be 
rather promising. 

(iv) The factorial moments can also be considered as a very sensitive sig- 
nature for clustering of particles in small bins of momentum phase-space. 
Indeed, a factorial moment of order q is sensitive only to the clusters con- 
tainig at least q particles. This obviously eliminates very effectively any back- 
ground. This point was recently illustrated by KLM collaboration analysing 
the collisions of 160 GeV/A Pb nucleons in emulsion [44]]. They found that 



events may differ drastically in the behaviour of their factorial moments, 
while no great difference is seen when they are analysed by other methods. 

(v) The distribution of the HBT radii obtained from individual events 
was also recognized since some time as a very interesting object to study. 
Recently, first data on this subject were presented by NA49 collaboration 
[]4"5| , [1[]. Although statistics is still limited (and the authors themselves do 
not attach too much meaning to the details of the plot) the results clearly 
show that the measurement is feasible and we may well hope for some exciting 
news in not-too-distant future. 



3 To study the moments of the factorial moments was suggested already some time ago 
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